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Phenomenological and numerical studies of the small scale spectra of energy are pre- 
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rH sented for high Reynolds number rotating Boussinesq flows in unit aspect-ratio domains 



^ We introduce a non-dimensional parameter r(A;) = jf^, where the wavevector k has 

rv^ vertical component fc^ and horizontal component fc/j, / is the Coriolis frequency and N 

> 

>0 is the Brunt- Vaisala frequency. While requiring that / and N are such that the potential 

cn 
in 



vorticity is nearly linear in the dynamical variables, we deduce that for F <C 1, the po- 
tential enstrophy suppresses the transfer of horizontal kinetic energy into large kh modes 
while forcing it to become independent of fc^, scaling as kj^^. For F 3> 1, the potential 
enstrophy suppresses the transfer of potential energy into the large k^ modes while forc- 
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^ ing it to become independent of kh, scaling as k^^. The usual scaling exponent —3 is 

adjusted down to —5 based on a heuristic argument which links the downscale fluxes of 
potential enstrophy and energy by their ratio. Spectra computed from high-resolution 
simulations of the Boussinesq equations with isotropic low-wavenumber forcing are used 
to explore such anisotropic constraints on the energy and provide a posteriori justifica- 
tion for the joint flux ansatz used to obtain the —5 scaling exponent. In simulations for 
which f/N <^ 1, the asymptotic regime of F ^ 1 is achieved and the horizontal kinetic 
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energy is shown to become independent of kz and scale as /c^'\ This collapse results 

in a uniform fc^ scaling for all k^ and hence a uniform fc*^ scaling of the energy in the 

spherical wavenumbcr, indicating extremely efficient net downscale transfer of energy. 

In simulations for which f /N 3> 1, the trend towards the asymptotic regime of F ^ 1 

is observed although the potential energy does not become entirely independent of fc/j 

as predicted. In simulations with f/N = 1 where rotation and stratification are equally 

strong, both F <C 1 and F ^ 1 regimes may be recovered depending on the aspect-ratio 

of the wavevector, t = kz/kh, which is also the cotangent of the wavevector angle to 

the vertical. The dependence of the spectra on r results, on average, in shallower scaling 

exponent values between —1 and —5/3 in both kh and kz implying an effectively more 

isotropic downscale distribution of energy in this limit. In all cases the empirical evidence 

points to both energy and potential enstrophy being jointly transferred downscale with 

the spectral scaling of the the former constrained by the latter. 



1. Introduction 

The Boussinesq approximation for hydrodynamic flow is the cornerstone for many 
applications in which a suitable model for scalar transport is needed. It is physically 
realistic in that it postulates that the fluctuations of the scalar are small compared to 
its background value and is computationally tenable for linear stratification profiles and 
simple boundary conditions. The Boussinesq approximation in the rotating frame forms 
the parent system for the derivation of further simplifications such as the shallow- water, 
hydrostatic and quasi-geostrophic (QG) models which are standard approximations in 
geophysical regimes which include rotation and stratification of the flow. 



High-wavenumber energy spectra of rotating Boussinesq 3 

The Boussinesq system for a stably stratified flow in a rotating frame is given by: 

—u + fzxu + Vp + N0z^vV^u + J' (f.f) 

— e-Nw^nW'^e (f.2) 

Dt ^ ' 

V-M==0, 

D d 

where u is the velocity, w is its vertical component, p is the effective pressure and T is an 
external input or force. The total density is given by pt{x) = po — bz + p{x), such that 
IpI ^ \bz\ <C Po where po is the constant background, b is constant and larger than 
zero for stable stratification in the vertical z-coordinate, p is the density fluctuation. 
The density in normalized to 9 ~ pig/bpo) ' which has the dimensions of velocity. The 
Coriolis parameter / = 20 where f2 is the constant rotation rate about the z-axis, the 
Brunt- Vaisala frequency N = {gb/po) ' , v = p/po is the kinematic viscosity and k 
is the mass diffusivity coefficient. We assume periodic or infinite boundary conditions 
with Prandtl number Pr = v/k ^ 0{1). The relevant non-dimensional parameters for 
this system are the Rossby number Ro = /„;// and the Froudc number Fr = fni/N, 
where /„; — (ejk'l)^'^ is the non-linear frequency given input rate of energy ej Smith & 
Waleffe (2002). Thus Ro and Fr are the ratios of rotation and stratification timescales 
respectively to the nonlinear timescale. 

There are three overlapping ideas from turbulence theory which we will exploit in 
this analysis of the Boussinesq equations. First from Kolmogorov (194f ) the notion of a 
universal range of scales governed by the dynamics of the flux of a conserved quantity. 
Second, the analysis of Kraichnan (f 971) on the inertial ranges of two-dimensional flows 
based on the constraint imposed by the joint conservation of both energy and enstrophy. 
Kraichnan thus deduced the k^'^ scaling, up to logarithmic correction, of the energy in 
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wavenumbers larger than the forcing wavenumber kf and k^^'^ scaUng for wavenumbers 

smaller than kf. Finally Charney (1971) used Kraichnan's ideas to show that similar 

constraints are placed on the quasi-geostrophic system due to the joint conservation of 

potential enstrophy and energy. Charney showed that in the quasi-geostrophic system, 

energy must flow upscale and potential enstrophy downscale, with predictable scaling 

exponents of the energy spectrum in the two ranges. 

The general principle then, is that for a given (inviscidly) conserved quantity in a 
system there is a range of scales over which the flux of that quantity is a constant. 
The larger the system and the more separation there is between the source (forcing) 
and the sink (dissipation, boundaries), the more extended this universal 'incrtial' range 
of scales becomes. Once this is established, it is fairly standard procedure to deduce 
phenomenologically what the distribution of such quantities must be in spectral space. 

In the inviscid and unforced {T = v = k = 0, ), Boussinesq system the potential 
vorticity g is a local (Lagrangian) invariant, as is its square the potential enstrophy Q, 
while the total energy is a global invariant. That is, 

potential vorticity q = (uJa- ^Pt ) , 7^ = 0, 

potential enstrophy Q = -q^ , = — / Q dx ^ 0, 

D f 

total energy Et = E + P, — Et dx ^ 0. 

At any point x, E — jlitp is the kinetic energy, P = ^9^ is the potential energy of 
the density fluctuations. The absolute vorticity Wa = a; + fz and the relative (or local) 
vorticity u> = V x it. PV may be written in terms of as 

BO 

q^ fN + io-V9+ f—^Nujs. (1.4) 

oz 

The constant part fN does not participate in the dynamics and we will therefore neglect 

it from now on. In what follows we will assume that j^ — > and k — > such that 
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Prandtl number Pr = v j k = 1, and the force T is confined to the lowest modes. Thus 

we assume a conventional 'inertial-range' of turbulent scales wherein the transfer of 

conserved quantities dominates over both their dissipation and forcing. Note that the 

global potential vorticity is zero for unbounded flows. While the potential enstrophy is 

also Lagrangian invariant, we will use its global conservation properties in our analysis. 

This set of invariants is analogous to those for 2-d turbulence where the vorticity is 

Lagrangian invariant and energy and enstrophy are globally conserved. It therefore seems 

reasonable to expect that an analysis similar to that of Kraichnan for 2D turbulence 

would yield similarly useful results. The raw Boussinesq system does not permit an exact 

relationship between energy and potential enstrophy as in the case of 2d turbulence. One 

notices that the potential enstrophy is, in general, a quartic quantity, while the energy is 

quadratic. However, when / and N are sufficiently large, PV linearizes in the dynamical 

variables and reduces in fourier space to: 

qik) ~ -i{fkj + Nkh X uh) = -lifkj + iNkhUh) (1.5) 

where ~ denotes fourier coefhcients, the total wavevector k = kh + k^z, the horizontal 
wavevector component has length kh ~ (fc^ + fc^)^'^, the vertical wavevector component 
has length k^ and it/j is the horizontal velocity vector with magnitude u^ — {u^ +u^)^'^. 
We assume that the vertical velocity w == u^ ^ in the leading order, consistent with 
empirical observations, and use incompressibility to obtain the last equality of Eq. (1.5). 
We will next show how in the formal limits where PV is linear, asymptotically exact 
relationships can be extracted between potential enstrophy and components of the total 
energy. 

First we brieffy review the theoretical results that have been obtained in the limits in 
which PV linearizes. For strong rotation and stratification, Embid & Majda (1998) and 
Babin et al. (2000) showed that the limit as Ro ~ Fr = e — t- is QG after the inertia 
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gravity waves are averaged out. The scalings predicted in Charney (1971) are recovered, 

namely k~^ scaling ol the energy for kf <$:. k <$:. kd, and k^^^^ scaling for k < kj where 
kf is the forcing wavcnunibcr and kd is the typical dissipation wavenumbcr. Embid & 
Majda (1998) also showed that the limiting dynamics for Fr -^ with finite Ro includes 
the so- called vertically sheared horizontal flows (VSHF). In both these cases, one can 
obtain a set of closed reduced equations for the linear PV dynamics from which all other 
quantities of interest can be derived. There is as yet no equivalent rigorous theoretical 
result for the unit-aspect ratio strongly rotating Ro — > case with finite _Fr, since in 
that limit the scalar equation 1.2, which is independent of Ro, vanishes in the leading 
order. This is a key reason why a (closed) reduced QG-like set of equations for this limit 
in the unit aspect-ratio domain has not been derived. Julicn et al. (2006) showed that 
if an additional small parameter is introduced, namely the aspect-ratio of the domain, 
then reduced equations are possible for the vanishing Ro limit as well. 

In recent theoretical work we predicted that the potential enstrophy associated with 
the Boussinesq rotating system exhibits universal statistical properties in the so-called 
'inertial range' of scales (Kurien et al. (2006)) in physical space for the same limiting 
regimes of the Rossby and Fronde mimbers described above. Formally, keeping lowest- 
order terms in a perturbation expansion in powers of e — > 0, universal statistics are 
found when the potential vorticity is linear in the dynamical variables u and 6. In those 
cases, the associated von Karman-Howarth equation von Karman & Howarth (1938) for 
the two-point correlation of q yields a simple scaling law for the third-order velocity- 
potential- vorticity correlation 

2 

{q{x)q{x + r){uL{x) - ul{x + r))) = -g^Q^ (1-6) 

where r is a physical separation scale, subscript L denotes the component along r, and 
Sq is the mean dissipation rate of potential enstrophy Q. The 2/3-law (1.6) assumes that 
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the two-point statistics of potential vorticity have a universal small-scale statistically 

isotropic component for statistically steady state flows. A related result was subsequently 

derived for the strict QG case (beginning with the evolution equation for qqa) Lindborg 

(2007) relating the third-order moments of potential vorticity difference and horizontal 

velocity difference across scale r to potential enstrophy flux, assuming axisymmctry of the 

small-scales, and restricted to components of the velocity that are in the horizontal plane 

(given that vertical velocity is strictly zero for pure QG). The numerical experiments to 

verify 1.6 are currently underway and will be presented in a separate work. 

There have been key numerical studies of the theoretical limits described above of which 

we here review the most relevant to our present study. The limit of strong rotation and 

stratification including the case f = N was examined by Bartello (1995) in a paper which 

studied geostrophic adjustment as a process in which wave- vortical- wave interactions 

acted as catalyst for highly efficient downscale transfer of energy downscale. Smith & 

Waleffe (2002) demonstrated the generation of slow large scales in strongly stratified 

turbulence with small scale forcing and the dominance of the fc/j = (VSHF) wave 

modes in the low wavenumbers. Waite & Bartello (2004) and Lindborg (2006) explored 

the effect of vortical forcing in stratified turbulence without rotation, thus highlighting 

the difference that 2D vs. 3D forcing can make in the evolution of the spectra. Sukhatme 

& L.M. (2008) showed that fairly small departures from f/N — 1 can have measurable 

effect on the scaling of the spectra for moderately low-resolution simulations. Waite & 

Bartello (2006a, 6) and demonstrated the transition to QG from stratified turbulence 

as Ro is decreased for fixed small Fr. Remmel et al. (2009) proposed a novel non- 

perturbative model to add non-QG effects to a flow and demonstrate the importance of 

resolving the small scales (waves) in order to obtain even gross large-scale features such 

as the energy growth accurately. 
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We will here study rotating Boussincsq turbulence both heuristically and using data 

from high-resolution simulations for Ro and Fr chosen such that the potential vorticity is 

nearly linear, and potential enstrophy is nearly quadratic. In our approach we retain the 

type of physically sound arguments used by Kolmogorov, Kraichnan and Charney, namely 

a focus on conserved quantities (in the inviscid, non-diffusive system), existence of inertial 

ranges where the downscale flux of a conserved quantity (in our case potential enstrophy) 

governs the dynamics and consequent spectral scaling of other related quantities (energy), 

and the notion of a statistically universal range of small scales. 

2. Spectral relationships bet^veen potential enstrophy and energy 

We begin by assuming that velocity fluctuations and density fluctuations are about 
the same order. We can construct spectral relationships between the quadratic potential 

enstrophy Q{k) = -|(7(fc)p and either the horizontal kinetic energy Eh{k) = -\uh{k)\^ 

1 ~ fk 

or the potential energy Pik) — —\0{k)\ based on the parameter T = — r^, which by 

2 Nkii 

definition resembles a local Burger number in spectra space. From Eq. 1.5, for F ^ 1, 

stratification dominates the PV q{k) ~ iNkhUh, (2.1) 

and potential enstrophy Q{k) — -\NkhUh\^^ (2.2) 

/•OO /'OO 

giving the constraint lim / Q{k)dkh > N'^h\ i Eh{k)dkh. (2.3) 

Eq. (2.3) is interpreted to mean that for wavevectors with large kh, that is the 'wide' 
wavevectors, the potential enstrophy Q{k) dominates the downscale (horizontal scales 
smaller than 1/ Kh) dynamics and consequently suppresses the transfer of horizontal ki- 
netic energy into those modes. This is analogous to the relation in 2d turbulence 

f-OO /'OO 

lim / Vt{k)dk :^ K^ I E{k)dk (2.4) 
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wherein the enstrophy Q{k) suppresses the downscale transfer of kinetic energy E{k) for 

sufficiently high wavenumbcr k, and in turn forces the inverse cascade of energy. 

In the opposite hmit F ^ 1 one can show that, 

rotation dominates the PV q{k) ~ —ifk^O, (2-5) 

and potential enstrophy Qik) — -;^\fkz9\'^, (2-6) 



/•OO /"OO 

giving the constraint lim / Q{k)dkz ^ /^^^ / P{k)dkz 



(2.7) 



which is interpreted to mean that potential enstrophy suppresses the transfer of potential 
energy P{k) into the wavevectors with large k^, that is the 'tall' wavemodes. In physical 
space this may be loosely interpreted to mean that potential energy transfer into the 
small vertical scales is vertical scales. 

More generally, for fixed f /N there is a transition regime when kh/kz — f /N if the 
wavenumber regime is large enough. Thus, in an idealized infinite wavenumbcr flow, for 
fixed f /N there will be some regime of wavevectors such that kh/kz ^ f/N for which 
r ^ 1 and horizontal kinetic energy is suppressed in kh according to Eq. 2.10, and a 
regime of wavevectors such that k^/kz ^ f/N for which F 3> 1 and potential energy is 
suppressed according to Eq. 2.11. The transition regime is defined by a cone in wavevector 
space with angle 6 to the vertical such that tan 6* = f/N. All wavenumbers inside the 
cone would behave according to F 3> 1 and all those outside would behave according to 
F <C 1. In our simulations in finite sized boxes, we study three extreme cases: f/N <^ 1 
where the cone is extremely narrow (small angle) so that most wavevectors lie outside 
the cone giving F <C 1, f/N ^ 1 where the cone is very broad (large angle) so that most 
wavevectors lie inside the cone and F ^ 1 and finally f/N = 1 where the cone-angle is 
45deg to the vertical and both F <C 1 and F :^ 1 regimes exist. 

There are several noteworthy differences between the relations Eq. (2.3) and (2.7) 
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derived above and the 2d relation Eq. (2.4). The latter relates two (inviscidly) conserved 

quantities in spectral space, while the former relates the spectrum of a conserved quantity 

{Q{k)) with that of one {Eh{k)) that is not conserved in general. The 2d turbulence 

constraint forces enstrophy downscalc and energy upscale. In rotating Boussinesq, for 

r ^ 1, when the horizontal kinetic energy is suppressed in the large kh, there is no 

corresponding constraint in fc^ for either the kinetic or potential energies; conversely for 

r 3> 1 which potential energy is suppressed in /c^, there is no corresponding constraint 

in kh- Therefore there is always the possibility of a downscale transfer of energy in the 

unconstrained wavevector component. That is, a dual cascade of potential enstrophy and 

energy is possible downscale. In addition, Eq. (2.4) is isotropic in wavenumber in that it 

depends on the total spherical wavenumber while our new relations Eqs. (2.3) and (2.7) 

are highly anisotropic in wavenumber. As we will show from the data, the spectra indeed 

indicate that energy is transferred downscale in all cases albeit in a highly anisotropic 

manner. 

2.1. A scaling phenomenology 

Following Kraichnan and Charney one can naively postulate a scaling exponent for the 
horizontal kinetic energy and potential energy. Assuming the potential enstrophy flux 
governs the downscale dynamics, the energy scaling in the high wavenumber regime 
must depend only on the potential enstrophy dissipation rate eq and the wavevector 
(component) magnitude. Dimensional analysis then results, for the two limits: 

r « 1: Eh{k) = CEe^^^'k^^ (2.8) 

r > 1: P{k) = Cpe'^'^k^^. (2.9) 

where Ce[p] is a prefactor which in principle should depend on the ratio of rotation to 
stratification since Eh and P are not separately conserved. In this study we focus on the 
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scaling in wavenumber and leave the study of the prefactor to a future work. Analysis 

of high-resolution data (see numerical results below and prior work in Kurien et al. 

(2008)) showed that the actual scaling exponent was much steeper than —3, prompting 

a refinement of the above scaling estimates as follows. 

We first suppose that part of the downscale energy transfer in both limits is not entirely 

suppressed. In the case F ^ 1 where we predict that Eh{k) suppressed for large kh^ it 

may be that Eh{k) or P{k) is free to transfer into modes with large k^- And in the case 

r :^- 1 with P{k) suppressed for large fc^, it may be that P{k) or Eh{k) is free to fill 

modes with large ku. Thus there may be a partial flux of the total energy (kinetic plus 

potential) to high horizontal or vertical wavenumber leading to the postulate that the 

spectral distribution of energy must depend on both eq and energy dissipation Sd in the 

small scales, where the subscript d denotes downscale energy. We will show a posteriori 

that our data supports such a hypothesis. A simple way to include the dependence on 

energy dissipation is to postulate that the energy scaling depends on the ratio Sq/ed and 

on the wavevector component. Then dimensional analysis yields: 

T^l:Eh{k)ozCE — k^^ (2.10) 

r>l: P(fc) ex Cp^/c-^ (2.11) 

£d 

Irrespectively of the exact value of the scaling exponent for the steep decay, for F <C 1 
the spectral horizontal kinetic energy distribution becomes independent of k^, while for 
F 3> 1 the spectral potential energy distribution becomes independent of kh ■ Notice also 
that this estimate requires that eg/e remains finite in the inviscid i^ — > 0, non-diffusive 
K -^ limit. One could presumably add corrections to mollify the blow-up for small fc/i 
and kz respectively along the lines of the log-correction derived by Kraichnan (1971). This 
issue of scaling correction will not be addressed in this paper except to say that steeper 
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than -3 scaling implies even stronger non-local transfers than 2d turbulence as discussed 

by Kraichnan (1971)). As we will show in the sections below, our modified assumptions 
about energy transfer downscale are supported by the results of our simulations, and the 
scaling exponents measured from our data for the appropriate regimes are remarkably 
close to —5. 

Bellet et al. (2004) studies mode-angle dependence of the total energy spectrum of 
rapidly rotating turbulence flow, using closure models for the anisotropic rotating equa- 
tions for wave turbulence. The flows they study do not have stratification or any advected 
scalar and thus have no potential vorticity. Therefore, our present analysis, which does 
in effect study spectra as function of mode-angle via /c^/fc/i, cannot be carried over to 
deduce similar constraints. What we will present poses several departures from Bellet 
et al. (2004). First, we address, in principle, a very broad parameter space of finite Ro 
and Fr at high Reynolds numbers via the parameter T. We also, for the first time, show 
that potential enstrophy conservation imposes constraints on the energy distribution in 
parameter regimes other than strict QG. And finally we propose a way to understand 
how a conserved quantity (potential enstrophy) may affect the spectral distribution of a 
non-conserved quantity (horizontal kinetic energy or potential energy) and furthermore 
propose a consistent argument to estimate the resulting spectral scaling exponents. 

3. Numerical simulations and results 

We perform pseudo-spectral calculations of the Boussinesq equations with rotation on 
grids of 640'^ points and 1024'^ points in unit aspect-ratio domains. The computational 
domain has length L = 1 to a side, with wavenumbcrs in integer multiples of 27r. The 
time-stepping is 4th-order Runge-Kutta with resolution of the smallest wave frequencies 
with five timesteps per wave period. In the simulations at 640"^, the diffusion of both 
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momentum and density (scalar) is modeled by hyperviscosity of laplacian to the 8th- 

power in order to extend the inertial scaling ranges. The Reynolds number in these 

cases is not computed explicitly but the hyperviscosity coefficient is chosen to resolve the 

total energy in the largest shell Chasnov (1994); Smith & Walcffe (2002) thus ensuring 

adequately turbulent ffow for a given resolution. The simulations at 1024"^ were run with 

regular Navier-Stokes viscosity. The energy input rate Sf ~ 0.5 and 1 respectively for the 

smaller and larger resolutions, and Ro and Fr are varied by varying the rotation and 

stratification rates for fixed energy input rate which determines the nonlinear timescale. 

The forcing is incompressible and equipartitioned between the three velocity components 

and 6. This is equivalent to forcing the two wavemodes and the one vortical modes equally, 

thus providing an optimally unbiased calculation. The forcing is also stochastic and the 

phases of each mode change sufficiently rapidly that the forced modes are statistically 

isotropic. The forcing spectrum is peaked at kf = 4 ± 1, for large scale forcing. The 

simulations are dealiased according to the isotropic two-thirds dealiasing rule. Lower 

resolution runs at 256'^ and 512^ (see Kurien et al. (2008)) were also performed in earlier 

studies, but are not reported here. 

We seek to verify our predictions for three particular "extreme" regimes in Ro, Fr 

(see Table 1). The first is strongly stratified with moderate rotation, f /N <C 1 (runs 1 

and 4); the second regime has strongly rotating with moderate stratification, f /N 3> 1 

(runs 2 and 5); and finally the regime which fixes f ~ N large, (run 3) so that the 

dependence of potential vorticity on the ratio F is reduced to a dependence on the 

cotangent t = k^/kh alone, of the angle of the wavevector with respect to the vertical. 

Presumably, for sufficiently long inertial ranges any combination of /, TV, fc/j and kz 

such that r 3> or <C 1, while maintaining near- linear PV, may be similarly investigated. 

As discussed in Bartello (1995), strict potential enstrophy conservation is not possible 
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in a truncated inviscid spectral scheme due to the lack of conservation of Lagrangian 

invariants by individual triads. The best we can expect with this method in the viscous 
case is that truncation errors are swamped by viscosity thus giving conserved fluxes in 
the inertial range. We mention this fact because we will be using the notion of potential 
enstrophy conservation in this study. The spectral method itself is routinely used to 
compute flows very similar to those we study here and and we proceed while remaining 
aware of the caveats. 

For each case in the Table, we will present the calculations to check specific predic- 
tions of the phenomenology proposed, namely the anisotropic suppression of horizontal 
kinetic energy and potential energy in the different F regimes, and the scaling exponents 
predicted in these regimes by Eqs. (2.10) and (2.11). The spectra as a function of {kh, k^) 
are defined as follows: 

l|a(fc)p = l|a(fc,,fc,)|2^i ^ \a{K,,k^r (3-1) 

fcft-0.5^Kh<fcfc+0.5 

where a is any of q, Uh or 9. We will present the horizontal kinetic energy and potential 
energy spectra (a) as a function of kh for various values of fixed kz , (b) as a function of 
fc/i, summed over kz, (c) as a function of kz for various values of fixed k^, and (d) as a 
function of kz, summed over fc/j. The (a)-(d) above correspond to subplots (a)-(d) in the 
sets of figures to follow. In each subplot, horizontal kinetic energy is plotted on top and 
potential energy on the bottom. We will also show the more conventional decomposition of 
the energy in contributions from the linear cigcnmodcs of the Boussinesq system. There 
are two inertial-gravity wave-eigenmodes with non-zero frequency and zero PV which 
contribute to the wave (or ageostrophic) energy spectrum. The remaining eigenmode has 
zero frequency and linear PV and contributes to the vortical (or geostrophic) energy 
spectrum. We will attempt to relate the mode-angle representation to the usual wave- 
vortical decomposition. Since this paper is largely an exploration of spectral scaling in 
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run grid diffusion e/ kf f (Ro) N (Fr) 

1 640* hyper 0.5 4 14 (0.48) 3000 (0.0023) 

2 640* hyper 0.5 4 3000 (0.0023) 14 (0.048) 

3 640* hyper 0.5 4 3000 (0.0023) 3000 (0.0023) 

4 1024* NS 1 4 8.58 (1) 8580 (0.001) 

5 1024* NS 14 4290 (0.002) 8.58 (1) 

Table 1. The parameters for the numerical simulations of the Boussinesq equations. The nota- 
tion "hyper" denotes hyperviscous diffusion of Laplacian to the eighth power, apphed to both 
momentum and scalar; "NS" denotes regular Navier-Stokes viscosity. 

the large wavenumbcrs, we will, for clarity's sake, use the hyperviscous simulations data 
with their extended ranges to compute scaling exponents. The Navier-Stokes cases (runs 
4,5) show no appreciable differences from their hyperviscous counterparts except for 
shortened scaling ranges, as is to be expected. Wherever appropriate, and when it does 
not result in a proliferation of figures, we will also present the data from the simulations 
with Navier-Stokes viscosity. 

First we show in Fig. 1 that in each of the three cases, the global potential cnstrophy 
(solid line) is captured by one half of the square of the corresponding linear PV, to within 
3% or better agreement, potential cnstrophy, after a period of spin-up and over-shoot, 
settles to a constant, indicating statistically steady state for this quantity, the mean rate 
of input of potential cnstrophy balancing the mean rate of dissipation. We are thus able 
to stablish the tinic-framc over which the small scales, governed by potential cnstrophy 
dynamics, have achieved statistically steady state. The small scale spectra of energy that 
we are interested in, achieve a statistically steady state. The total energy meanwhile (not 
shown) continues to grow since it can transfer upscale, where in our case, there is no sink 
to drain the energy. The calculation of all statistical quantities in what follows occurs 
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Figure 1. [Color online] Global potential enstrophy evolution in time for runs 1,2, and 3. The 
total Q and its quadratic components corresponding to linear q, are shown. Solid line: total 
Q\ dashed line : ll/fj + Nijj-j,\^\ stars: ^\N ijj-j,\^ \ circles: ||/fjP- In each limiting case, the 
corresponding quadratic part of the potential enstrophy captures the total Q to within 3%. 

over the time-frame where the m.ean potential enstrophy is constant, 3.5 < t < 6, over 
about 25 frames. 



3.1. Strong stratification, moderate rotation 

This case is studied using data from runs 1 and 4 of Table 1. We are in the limit F ^ 1 
since nearly all wavevectors lie outside the narrow cone with vertical angle tanO — f /N = 
0.0047. Therefore, q{k) = q{k}i,kz) ^ Nk^Uh irrespective of k^/kh and the spectral 
scaling of Eh{k) should become independent of kz as given by Eq. 2.10. Fig. (2(a)) 
(top) and shows that Eh{kh,kz) collapses to a function of fc/i alone over the wide range 
< A:^ < 100, and scales very close to fc,^^. In fact, the collapse is so complete that the 
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(a) Spectra of horizontal kinetic energy 
(solid lines) and potential energy (dashed 
lines) as a function of kh for various fixed 
1 ^ fch ^ 100 as indicated. 



(b) Spectra of horizontal kinetic energy 
(solid line) and potential energy (dashed 
line), as functions of kh, after summing over 
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1 ^ fch ^ 30 as indicated. kz- 

Figure 2. [Color online] The horizontal kinetic energy and potential energy spectra for f <^ N 
moderately rotating and strongly stratified flow (Table 1 run 1). 
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converse plotting of Eh as a function of k^ in Fig. (2(c)) for various fixed kh, gives a series 

of plateaus, indicating no k^ dependence at all in the inertial range for any particular 
kfi- The corresponding plots for P{kh,kz) shows that it nominally follows the scaling 
of Eh as a function of k^ (Figs. 2(a)), yielding /c^-averaged scaling of fc^^ (Fig- 2(b))). 
However the converse plotting of P (Fig. 2(c)) shows that it does not become completely 
independent of k^ even though, when summed over kz (Fig. 2(d)), its scaling becomes 
close to /c^'' as predicted in this limit for Eh{kh). Thus P appears to mimic the scaling 
of Eh as a function of kh, but not as a function of kz- 

We now make a more conventional representation of the spectra by plotting the en- 
ergy contained separately in the wave and vortical modes as a function of the spherical 
wavenumbcr k. Figure 3 shows that almost the entire downscale energy is captured by the 
vortical modes. The narrow upscale range is also filled by the vortical modes. The total 
energy (vortical energy) distribution is nearly a constant as a function of k. The wave 
energy is subdominant everywhere except in the forced modes, indicating that the wave 
energy does not move efficiently out of those modes in this strongly stratified regime. This 
figure is to be compared with Fig. 2(d) to see that for k > kf the vortical modes mimic 
the horizontal kinetic energy and that the independence of the latter from kz translates 
to a constant spectral distribution over k as well, given that the dependence of Eh on kh 
is rapidly very decaying (Fig. 2(b)). To the best of our knowledge, this is the first obser- 
vation and explanation of the asymptotic high wavenumbcr energy spectrum for strongly 
stratified flow with flnite rotation, with isotropic forcing in the low-wavenumbers (large 
scales). Waite & BartcUo (2004) observed fc^^ scaling for the total energy spectrum and 
the vortical mode energy spectrum in flows with rather low stratiflcation N, zero rotation 
and forced in the vortical modes only. They also observed saturation in kz of the total 
energy spectrum. In the light of the analysis in the present work, it would appear they 
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Figure 3. [Color online] Strongly stratified and moderately rotating cases, with hyperviscosity 
(run 1, left) and with Navier-Stokes viscosity (run 4, right). Total (black solid), wave (red dash), 
vortical (blue dot-dash) and fch = mode (blue dotted) energy spectra. Almost the entire total 
energy is captured by the vortical modes within which the energy is uniformly distributed as 
fc". The wave-mode contribution dominates in the forced modes only. The fch = modes form a 
negligible contribution, highlighting the difference from VSHF. 

achieve something like the regime F <C 1 and linear q ~ A^cjs, not with large A'^, but with 
large ujj, achieved by forcing vortical modes only. Their paper Waite & Bartello (2006 &) 
which studies similar parameter regimes in Fr while varying Ro does show fc" scaling of 
the vortical mode energy for Ro ^ 0(1) but does not recover the k^^ steep decay. As 
we have pointed out above, the exact decay exponent does not matter as long as it is a 
decay and is independent of k^ , leading to the fc^ scaling. Our analysis goes a step further 
in proposing an explanation of why there is suppression in fc/j and and not in A:^ due to 
the very particular constraints imposed on energy by the potential enstrophy, leading to 
independence from kz- We have also conclusively shown that the asymptotic small Fr 
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regime for the small scales can be achieved with unbiased isotropic forcing and in the 

presence of finite rotation. 

There is also a stark difference from the vertically sheared horizontal flows (VSHF) 

{kfi = modes) expected to be included in this limit when the forcing is at high wavenum- 

bers (small scales) . In our data the kh = modes contain a very small percentage of the 

total energy (see Fig. 3) and the vortical (PV) modes contain most of the energy, whereas 

the opposite is true for the low wavenumbers in the case where the forcing is in the high 

wavenumbers (Smith & Waleffe (2002); Embid & Majda (1998)). 

3.2. Strong rotation, moderate stratification, Ro = 0.0022, Fr = 0.48 

For run 2 of our simulations with moderate stratification and strong rotation, F 3> 1 and 
we expect to observe Eq. (2.11). In Fig. 4(a) we see that although the scaling of P(k) in 
kz for all kh shows some degree of agreement, the curves themselves have not collapsed. 
This is also shown by the converse plotting of P{k) for various k^ functions of k^ in 
Fig. 4(a). However we note that the trend as Ro decreases is for the potential energy 
curves to collapse to functions of kz alone. Thus Ro — 0.002 may not be sufficiently 
small to see the limiting small-scale dynamics for Ro — > 0, Fr ^ 0{1). And this despite 
the fact that potential enstrophy has converged to the stratification dependent quadratic 
component as seen in Fig. 1. Thus it appears, at least empirically, that for the small 
scales, the asymptotic regime of strong stratification is easier to achieve than that of 
strong rotation. This would seem related to the fact that the Ro -^ limit for finite 
Fr is not amenable to a reduced set of equations since in that limit the scalar evolution 
equation vanishes in the leading order. 

Figures 4(b) and 4(d) show that the fc^-summed spectra scale between k^ and fc^ 
and the fc/i-summcd quantities scale approximately as k~^. Thus, on average both hori- 
zontal kinetic and potential energies are suppressed in the large kz modes while both are 
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free to fill the large kh modes. Again this justifies the ansatz that the scaling depends 

on eq/ed, with downscale fluxes of both potential enstrophy and energy fluxes playing a 

role in the final development of the high wavenumber spectra. 

The conventional representation of wave and vortical mode spectra, Fig. 5 shows that 
in this regime, the vortical modes once again dominate the contribution to the total en- 
ergy at all wavenumbers except around the forced wavenumbers, where the wave- mode 
contribution dominates. At high wavenumbers, the vortical modes scale more or less as 
/j-s/s ^}jgQ compared with Fig. 4(d) this may be interpreted as the vortical modes being 
dominated by the potential energy, which persists strongly into the large k^ modes while 
decaying rapidly in the large kz modes. Thus in the aggregate the k^ scaling dominates in 
the spherically averaged spectra as well. On the right of Fig. 5 shows the same decomposi- 
tion for the run 5 in the strongly rotating regime which shows qualitative agreement with 
the hyperviscous case, comparable scalings and hence providing some confidence that in 
these cases hyperviscosity does not bias our interpretation with undesirable artifacts. 

The point of similarity between the strongly stratified versus the strongly rotating 
limits is that in both cases the vortical modes dominate the high wavenumbers or small 
scales. However, in the former the vortical modes carry the horizontal kinetic energy 
whereas in the latter the vortical modes carry the potential energy. This is very different 
from QG wherein the vortical mode energy decays as k~^ in the high wavenumbers. 
Note also the difference from the equally strongly rotating and stratified case (see below) 
which, as is well known (see Bartcllo (1995)), carries most of the high wavenumber energy 
in the wave modes which scale as k~^. 

3.3. Equally strong rotation and stratification, Ro = Fr ^ 0.0022 

In the special case oi f ~ N the dependence of PV on F reduces to a dependence on the 
cotangent r = k^/kh of the wavevector angle. As r <C 1, the TV-dependent part prevails 
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while for T 3> 1, the /-dependent part takes over. That is, the cone in fc-space separating 

r <C 1 and r ;» 1 makes a 45° angle to the vertical allowing equal numbers of modes 

in both regimes. Using run 3, we show that as t <C 1 the value of \q{k)\ ~ N\khUfi{k)\ 

(Fig. 6, left panel), while as t ::^ 1 the value of |<7(fc)| ~ f\kzd{k)\ (Fig. 6, right panel). 

Since there is thus an average over the annulus of horizontal wavcnumbcrs contained in 

kfi ± 0.5 the curves on the left panel of Fig. 6 are smoother than those on the right. The 

special case of fc/j = has been omitted for clarity on the right panel as it is noisy (less 

averaging), while the k^ ~ case on the left poses no such problems. This figure clearly 

shows the dependence of q{k) ~ q{kh, k^) on the mode-angle with respect to the vertical 

as expected by Eq. (1.5) for f = N large. 

Next we consider the effect on the scaling of energy spectra as predicted in Eqs. 
(2.10,2.11). Fig. 7(a) shows E/^ and P as a function of kf^ for various fixed k^- For Ei^, the 
spectrum for a particular k^ changes over from basically flat to slowly increasing for small 
kfi to close to k'f^' decay when kh > kz- This is in agreement with the deduction that 
horizontal kinetic energy is suppressed in the small aspect ratio modes by the dominance 
of potential enstrophy in those modes, and consequently scales according to Eq.(2.10). 
Note that there is a cross-over range below which the energy decreases as k^ increases, 
and above which the energy increases as k^ increases. In Fig. 7(a) the spectra for po- 
tential energy for which we do not have a prediction in this limit of aspect ratio, when 
plotted in the same manner as above, also shows aspect-ratio dependence for r <C 1 but 
with much shallower scaling, the fcj^ line is shown for comparison. There is no cross-over 
range for P{kh, k^) and there appears to be much less suppression of potential energy for 
modes such that r <C 1. 

In Fig. 7(b) we show Eh{kh) summed over all kz (top) and similarly P{kh) (bottom). 
Interestingly, these scale approximately as fc^ and fc^ respectively, which is much 
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shallower than the fc^' observed or Eh for r <C 1, and even shallower than fc^ . Thus, 

while for a given k^ the horizontal energy is suppressed for r <C 1, the average over kz 

indicates net transfer to large fc/i of both Eh and P. 

To verify the prediction Eq. (2.11) for the large aspect-ratio modes r ^ 1 we examine 
the scaling of P{k) as a function of kz for various fixed kh (Fig. 7(c)). Here again the 
spectrum is nearly a plateau with a turnover to /cj^ as t ^ 1. Note again the cross- 
over range below which potential energy decreases as kh increases, and above which the 
potential energy increases as kh increases. Eh{k) plotted in the same way also shows 
some aspect ratio dependence but with much shallower scaling indicating that Eh is not 
suppressed in the wavevectors with t 3> 1. 

In Kurien et al. (2008) we presented the spectral scaling results for the case f = N 
for runs at 512^ and all other parameters very similar, except that in that case only 
the momentum was forced. In that paper, we did not have an estimate for steeper than 
—3 scaling exponent. The resulting spectra behaved very similarly to what we have 
presented in Fig. 7 so far as scaling exponents are considered. The overall magnitude of 
the potential energy spectra are much higher in the scalar-forced case presented here, as 
is to be expected. Thus the detailed forcing seems not to affect the results as far as the 
scaling exponents go. 

To conclude this discussion for the f ^ N strongly rotating and stratified case, the 
suppression of Eh in the small aspect ratio modes k with r ^ 1 does not result in 
potential energy being suppressed in the same way. Conversely, the suppression of P in 
the large aspect ratio modes with r 3> 1, does not result in the horizontal kinetic energy 
being suppressed in those modes. However, on summing over all kz [kh] because of the 
strong aspect-ratio dependence, Eh{kh) [P{kzy\ scales with exponent between —1 and 
—4/3. Thus a mode-by-mode suppression of horizontal energy in the small aspect-ratio 
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modes and of the potential energy in the large aspect-ratio modes, nevertheless results 

in a net transfer of these energies into the high wavenumbers on average, along with the 

downscale transfer of potential enstrophy. We thus obtain a posteriori justification in our 

ansatz of eg/ed dependence (Eqs. (2.10) and (2.11)) for the functional form of E^ and 

P. 

This example also goes to show how one-dimensional spectra in either the k^ or kz 

wave-component may obscure underlying dependencies on aspect-ratio of the modes and 

hide the subtle ways by which the energy populates the high wavenumbers on average. 

This possibility was noted in Bellet et al. (2004) who studied purely rotating flow, with 

no PV, and noticed a dependence of the spectra on the mode-angle as well. 

4. Flow visualizations 

To give a qualitative picture of the quantities of interest, we present volume rendering 
of snapshots of the fully developed flow for the horizontal kinetic energy and the potential 
energy in Fig. 4. It must be noted that it is not clear what correspondence the spectra have 
to the structures in the physical space visualizations of the flow. The spectra are related 
to the physical space correlations as fourier transform pairs. A particular distribution 
of a field in physical space does not indicate a unique spectral distribution in fourier. 
The visualizations presented here should be treated as qualitative impressions of the fiow 
structures. Nevertheless, we will point out consistencies between what is observed in the 
spectra with what is observed in visualizing the corresponding physical space quantities. 

When stratification dominates over rotation the horizontal kinetic energy exhibits flat 
horizontal structures at all values (Fig. 9(a)) which appear to have very little extent in 
the vertical direction. This is consistent with both Eq. 2.10 and the discussion and plots 
in section 3.1 which confirm complete independence from k^. These structure are quite 
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different from VSHF since the kh = modes have a very small contribution in this regime. 

The potential energy for this flow (Fig. 9(b)) appears to exhibit some tendency towards 

horizontally oriented blobs but retains more three-dimensionally coherent structures at 

the larger values, indicating that the potential energy distribution follows the distribution 

of horizontal kinetic energy to a degree, consistent with Fig. 2(b) but still retains some 

dependence on kz for small k^ (large vertical scales) consistent with Fig. (2(d)). 

Conversely, when rotation dominates over stratification. Fig. 9(d) shows that the po- 
tential energy exhibits structures that are strongly oriented in the vertical direction at 
all values with very little extent in the horizontal, consistent with Eq. (2.11). The hor- 
izontal kinetic energy also shows vertically coherent structures Fig. 9(c) but also shows 
high-valued larger blobs with more horizontal extent showing that E^ also retains some 
horizontal scale dependence. 

In the case of equal rotation and stratification (Figs. 9(e), 9(f)), horizontal kinetic 
energy and potential energy look rather isotropically distributed with no preferred direc- 
tion, in agreement with the spectra in Figs. 7(b) and 7(d) which show relatively shallow 
(ranging between -1 and -4/3) scaling in both k^ and kh on average. 

To summarize, equal rotation and stratification docs not show any preferential direction 
for the energy isosurfaces, consistent with nearly the same shallow scaling in k^ and 
kz (Figs. 7(b), 7(d)). Stratification dominated flows align the energy isosurfaces in the 
horizontal direction while rotation dominated flow align the energy isosurfaces along 
the vertical. The strongly stratified data are more in agreement with our asymptotic 
predictions than the strongly rotating data; the latter observation could indicate that 
even stronger rotation is required to see the asymptotic behavior. 
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5. Discussion and Conclusions 

The small scale, high wavenumber spectra of rotating Boussinesq have never before 
been studied at such asymptotically low values of Ro and Fr using high-resolution sim- 
ulations with unbiased isotropic forcing and dissipation. The rigorous theoretical results 
for the Ro ^ 0(1), Fr -^ by Embid & Majda (1998) and the corresponding simulations 
of Smith & Waleffe (2002) extract the kh = wave-modes as the VSHF contribution to 
the large scales. In our study of the small scales in this limit, which we have checked have 
negligible contribution of the fc^ = modes (see Fig. 3), the vortical modes display a 
constant spectrum for k > kf, apparently dominated by the fc" behavior of the horizontal 
kinetic energy (see Fig. 2(d)). This appears to be a consequence of one of the key new 
results of this work, namely that the horizontal kinetic energy becomes independent of 
(constant in) k^ for any kh (see Eq. 2.10 and Fig. 2(c)) for F <C 1. This uniform spectral 
distribution implies a very efficient downscale transfer, with energy contained predom- 
inantly by the vortical modes. The exact value of the steep scaling exponent in k^ is 
not as important for this result as the fact that the spectrum becomes independent of 
(constant in) k^ while decaying sufficiently quickly in k^ such that the spectrum becomes 
independent of spherical wavenumber k as well (Fig. 3). 

The strongly rotating case with finite stratification (runs 2 and 5) were expected to 
achieve the limiting behavior of Eq. 2.11 with collapse of the spectral curves for all k^ as 
a function of k^. While we see the correct trend in this direction, the curves in Fig. 4(a) 
show that the collapse is not complete. We are thus lead to believe that Ro will need to 
be driven even lower and the resolution increased to see the collapse. We can postulate 
that, for sufficiently small Ro, the potential energy curves for various kh will collapse to 
the same function of fc^, corresponding to Eq. 2.11. Then, potential energy as a function 
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of k would be flat, governed by the flat spectrum in kh and strongly decaying spectrum 

in kz- Again the exact exponent of the decay law would not affect this general conclusion. 

The equally strongly rotating and stratifled case (run 3) show how underlying anisotropic 

scalings of the spectra, depending on the cotangent of the mode-angle r, can be masked 

by summing over one or other of the wavevector components, as is commonly done. In 

Fig. 7(b) we show that the sum over all kz of the horizontal kinetic energy Eh{kh) shallow 

4/3 

with scaling ^ fc^ . Thus the average over all kz smears out the effect of suppression 
by potential enstrophy and results in a much shallower scaling that the r dependent 
fc^ scaling. The net effect is that energy flUs the high-wavenumbers with almost the 
same shallow scaling in kh and kz (Fig. 7(b) and 7(d)). In Fig. 8 the wave- vortical mode 
decomposition of the spectra shows that the total downscale energy scales as k^^ and is 
dominated by wave modes. This is a well-known result and the scaling is identical to that 
of the passive scalar since the wave-modes are passively advected in this limit. What we 
further establish from this analysis is that the shallow spectrum is set up in spite of the 
underlying steeper, anisotropic, r-dependent spectra in kh and kz- 

We have made an attempt at a physically consistent phenomenology that results in a 
steeper than -3 scaling of the spectra within this framework. The —5 scaling that arises 
out of the assumption that the ratio SQ/ed, where Cd is the net flux of the energy down- 
scale, plays a role in the downscale dynamics seems justified a posteriori. The analysis 
and verification of such an ansatz requires a detailed study of the fluxes and dissipation 
of both potential enstrophy and energy and is the subject of a separate work. As we have 
said before, the main conclusions of this paper hold irrespective of the speciflc scaling 
exponent for the decay of energy in the k^ and kz ■ 

The prediction for steep decay in kh according to Eq. (2.10) arises from the constraint 
that potential enstrophy places on the horizontal kinetic energy in kh for F <C 1. While 
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there is no explicit constraint on the potential energy, the data analysis in this regime 

reveals that potential energy scales very similarly to the horizontal kinetic energy as a 

function of kh, indicating that it too is suppressed in kh- Thus the total energy, where the 

contribution from the vertical component of the velocity is measured to be very small, 

appears to be suppressed in the large k^ modes. The observed shallow scaling in k^ of 

both horizontal kinetic and potential energies, suggests that the total energy is allowed to 

transfer into the large A:^ modes. Similarly, for F ^ 1, the suppression of potential energy 

in fcz is mirrored to a large degree by the horizontal kinetic energy. These observations 

indicate that the exchange between kinetic and potential energies is local in wavenumber 

in both f/N <C 1 and f/N ^ 1 regimes. 

Recent work explored the connection between hyperviscosity, Galcrkin truncation and 
the energy bottleneck Frisch et al. (2008) in Navicr-Stokes turbulence without rotation or 
stratification. In that paper it was proposed that high values of hyperviscosity result in 
a partial thcrmalization of the high wavenumbers which manifests as energy 'bottleneck' 
phenomenon which contaminates the inertial range. In the present studies of rotating 
and stratified turbulence, the only substantial difference between the hyperviscous and 
the comparable runs with regular NS viscosity is the consistently less extended scaling 
ranges of the latter runs, which is to be expected. Therefore we use the hyperviscous 
runs to extract scaling exponents while checking against the NS runs to make sure that 
bottleneck-type artifacts are not introduced. It appears that for rotating and stratified 
flow and for the value of hyperviscous power used, the bottleneck-type artifacts, if any 
arc minimal, and do not affect the purposes of this study. 

While we have discussed three extreme parameter regimes, this discussion encompasses 
all values of F thus allowing us to analyze arbitrary f/N in the limit of linear (or near 
linear) PV. In an idealized infinite domain where k — > oo, for any fixed J /N the F ^ 1 
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regime is achieved when kh/k^ 3> f /N while the F ^ 1 regime is achieved when k^/kz ^ 

f /N . The transition between regimes occurs for F = 1 or when 1/t — k^/kz ~ f/N. 

Thus we observe that this type of analysis could be very relevant to geophysical flow 

regimes where / and N values are not as extreme but where the scaling ranges are much 

longer. 
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Figure 4. [Color online] The horizontal kinetic and potential energy spectra for f ^ N strongly 
rotating and moderately stratified flow (Table 1 (ii)). The quantities plotted for subfigures 
(a,b,c,d) are the same as those for Fig. 2(c,d,a,b) respectively. 
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Figure 5. [Color online] Strongly rotating and moderately stratified cases, with hyperviscosity 
(run 2, left) and with Navier-Stokes viscosity (run 5, right). Total (black solid), wave (red dash), 
vortical (blue dot-dash) energy spectra. Almost the entire total energy is captured by the vortical 
modes within which the energy is uniformly distributed approximately as fc^^'^. The wave-mode 
contribution dominates in the forced modes only. 
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Figure 6. [Color online] Run 2,, f = N = 3000 {Ro = Fr = 0.0022). Left: \q{k)\ normalized 
by \khUh{k)\ as a function of kh for various values of kz- The normalized potential vorticity 
asymptotes to A'^ = 3000 for V = t <^ 1 as predicted in Eq. (2.1). Right: |g(fc)| normalized 
by \kz9{k)\ as a function of fcz for various values of kh- The normalized potential enstrophy 
asymptotes to f — 3000 as F = r 3> 1 as predicted in Eq. (2.5) 



High-wavenumber energy spectra of rotating Boussinesq 



35 



10* 


^^=^ 




... f|^'y 

k .30 


xio" - 


10' 


[ZIII^^^^^^ 




- 




zzi^^^^^^^~^^^^^^ 


10" 


iC^ 


'^ 


^5?^ 


' 


10"" 
10'' 






X^ 


% 




20 ,-''■' _---'" 


k- 


s555»»^ 




10 
-a 


30- -- 






"^ 









— EA) 


10"' 


Cv_^ 




...P{k„l 


10"' 


~~~\^ 


c 


- 


ID"" 


< 


'^ 


•-^ 


10"* 






\ '■• 



(a) Spectra of horizontal kinetic energy 
(solid lines) and potential energy (dashed 
lines) as a function of kh for various fixed 
1 ^ fcz ^ 30 as indicated. 



(b) Spectra of horizontal kinetic energy 
(solid line) and potential energy (dashed 
line), as functions of kh, each summed over 
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(c) Spectra of horizontal kinetic energy (d) Spectra of horizontal kinetic energy 

(solid lines) and potential energy (dashed (solid line) and potential energy P{kz) 

lines) as a function of kz for various fixed (dashed line), each summed over fch, as func- 

1 ^ fch ^ 30 as indicated. tions of kz . 

Figure 7. [Color online] The horizontal kinetic and potential energy spectra for f = N 
strongly rotating and stratified flow (Table 1 data (i)). 
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Figure 8. [Color online] Wave, vortical and total spectra for run 1, Ro — Fr = 0.0023. 
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